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Abstract 

Contrary to some previous claims, we find a sizable instanton contribution to the finite 
energy sum rule used to extract the value of the strong coupling from the measured 
r decay widths. It is of the same order of magnitude as standard nonperturbative 
corrections induced by vacuum quark and gluon condensates. Our result indicates that 
there might be no hierarchy of power corrections in finite energy sum rules at the 
scale of t mass. Therefore, the standard nonperturbative corrections do not necessarily 
improve the accuracy of the theoretical predicition, but can rather be used to estimate 
an intrinsic accuracy of the pure perturbative calculation, which turns out to be rather 
high on this evidence, of order one percent. 



*On leave of absence from St. Petersburg Nuclear Physics Institute, 188350 Gatchina, Russia 



1. In recent years the hadronic decay widths of the r lepton have attracted considerable 
attention []lj, p[ ^j] . The main interest originates from the observation that the inclusive 
nature of the decay widths allows for a theoretical prediction, comparison of which with 
the data can fix the value of the strong coupling with an accuracy competitive to its 
determination at the Z° pole. A direct measurement of the coupling at the r scale 
might provide most conclusive evidence for the perturbative running of the effective 
QCD coupling. This is important, since the values of the coupling extracted from deep 
inelastic scattering data and at the Z° pole could indicate a slower evolution compared 
to the QCD renormalization group equations, a possibility that has already triggered 
some activity towards alternatives to the standard QCD evolution |5], || [7|. 

Due to momentum analyticity, all dynamical information on r decays is contained in 
the two-point correlation functions of flavour- changing vector and axial vector currents, 
which are analyzed in the framework of the SVZ operator product expansion ||. In 
this approach the perturbative expansion of the correlation functions is complemented 
by nonperturbative corrections, which are proportional to the vacuum expectation val- 
ues of local operators build of quark and gluon fields (condensates) and which are 
suppressed by increasing powers of the r mass m T ~ 1.8 GeV. The main evidence for 
the applicability of standard SVZ expansions comes from QCD sum rules, which are 
believed to work with the typical accuracy of 10-20%. The problem is that a determi- 
nation of a s {m T ) with, say, 10% accuracy requires knowledge of the decay widths at 
the one percent level. It has been found that "standard" nonperturbative corrections 
to the r decay rate do not exceed a few percent Q. This still leaves open the possibility 
of nonperturbative effects that are beyond the SVZ expansion, an issue that should be 
investigated separately. 

In the recent paper || Nason and Porrati study the contribution to the total r 
hadronic width from small size instantons in the QCD vacuum. They find that even 
with the highest value of a s (m T ) allowed, the instanton contributions never grow beyond 
10~ 6 — 10~ 5 , and, therefore, are completely negligible. This smallness can essentially 
be traced to the well known fact ]10| that the instanton density is proportional to the 



product of light quark masses m u mdm s and instanton-induced transitions vanish for 
strictly massless quarks. This is only true, however, for one-instanton contributions in 
an empty (perturbative) vacuum, which is far from a physical reality. As explained 
at length in [ll], the density of small instantons of size p in the background field of 
large-scale vacuum fluctuations is modified in such a way that the current quark mass 
is substituted by an effective mass, 



This effect of the vacuum "medium" is very large. To obtain a rough estimate, one may 
replace the product of current quark masses in the answer given in |J by the product 
of constituent masses M u ~ ~ 350 MeV, M s ~ 500 MeV. This yields an increase 
by a factor ~ 10 4 , and boosts the contribution obtained in || into the range of percent 
corrections, which are important. Thus, a more quantitative analysis of this effect is 
mandatory, and this is the subject of the present paper. More precisely, we calculate the 
instanton-induced contribution to the relevant two-point correlation functions, which is 
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related in the framework of the operator product expansion to an exponential correction 
to the coefficient function in front of six-quark operators {qq) 3 - Since the current quark 
masses on the right hand side of ([I]) are small compared to the second term, we shall 
omit them altogether in the following. 

2. The correlation function of interest is the T-product of two flavour-changing vector 
or axial vector currents, 

U;i(q) = i J dAe^ A <0|r(j^)>(y)) |0) A^x-y, (2) 
j^x) = u(x)j^ \ \ d(x) , (3) 



of two quark fields, which we denote by u and d, in a theory with n/ massless quarks. 
At large euclidean q 2 the characteristic size p of the instanton scales as p 2 ~ q~ 2 <C A -2 , 
where A is a typical QCD scale parameter. Therefore an effective Lagrangian approach 
is appropriate. To lowest order in (pA) 2 an instanton of size p (in the singular gauge) 
induces an effective 2n /-quark vertex with the quark legs corresponding to the instanton 
zero modes. Since these zero modes possess definite chirality it is most convenient to 
use the two-component Weyl spinor notation: the euclidean quark fields and 7-matrices 
are written as 

9 = ( f a ) . * = > > = ( a l '7 ) , (4) 

where k = 1, 2, . . . , N c is a colour and a, a = 1,2 are spinor indices. We use the 
notations a^ a = (1, — ia), cr^a = (l,ia), where a are the Pauli matrices, and for 
vectors define v = v^o^v = v^a^. Whenever a mixture of spinor and colour indices 
takes place, it is understood that spinor matrices act in the 2x2 upper left corner of 
N c x N c colour matrices. The instanton vertex is then described by 't Hooft's effective 
Lagrangian 

n f 

L*„ = {40*3)^01, Oi = H(xi<p)(RiPi), (5) 

i=i 

where <p ka = e al3 Ujp, Rk a = ~ e ap{U\)ki Uj is the ^[/(iVcJ-matrix of the instanton 
orientation and e a/3 is the antisymmetric tensor with e 12 = 1. Colour and spinor indices 
will be suppressed, whenever no confusion can arise. 

Let us recall the derivation of this expression. To find the coefficient multiplying the 
operator O/, one should compare the 2n/-quark Green function IT^i i , i{ x i)Xi{Vi)i evalu- 
ated in the instanton background in the near mass shell limit Xi, yi — * 00, with the result 
obtained from the effective Lagrangian (|5|). The 2n/-quark amplitude is simply obtained 
by substituting the zero modes for the quark fields, i.e. equals Jliii K o(a ; fc)^o(yA,0- Ex- 
plicit expressions for the zero modes for the instanton with center at the origin are 
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xip 2irp 3 / 2 _ kx 2irp 3 / 2 

where = 1 + p^/a; 2 . For the instanton with center at Xo, the zero modes are obtained 
by the obvious substitution x — > x — x . The large distance behaviour of the zero modes 
coincides (in the singular gauge) with that of the perturbative propagator up to the 
factors 27rp 3 / 2 (p, 27rp 3 / 2 K, respectively. Therefore the leading term of an expansion of 
the zero modes in p 2 /x 2 is indeed reproduced by an insertion of the effective Lagrangian 
(H). Subsequent terms in this expansion correspond to the exchange of soft gluons, that 
is, to higher dimensional operators in the effective Lagrangian containing extra gluon 
fields. 

Let us now return to the correlation functions. In two-component notation 

j^x) = Mx)a^ d {x) + \Xu{x)<y^Xd{x) , A = j ~\ ^vector ' ( 7 ) 

In order to find the one-instanton contribution to the coefficient function of the operator 
of lowest dimension, Oj, one inserts the current product into a 2n/-point Green function, 
evaluates it in the instanton background and takes the large distance asymptotics. The 
relevant diagrams for all possible contractions of 



n ^( z *)x(4) 



Li=l 



jU x )jv(y) zi,4^°o (8) 



are shown in Fig.l, where the instanton is depicted as its effective 2nj-quark vertex. 
The solid line denotes the quark zero modes and the dashed line represents the quark 
propagator over the nonzero modes in the instanton background. Its explicit form is 
known from |T3|: 



2 U I xyJj\\ a^pj UixAa^yU] \ 

l + Pl 12^2 + 



Si(x,y) = i(x(x)x(y))i 




x 2 y 2 I Air 2 A 2 x 2 y 2 (y 2 + p 2 ) 



A / 2 U lX yU\\ , v»P 2 U lX a^AyU\ 

+ o_2 A 4 1 + Pi „2..2 + 



2vr 2 A 4 I 11 x 2 y 2 J Air 2 A 2 x 2 y 2 (x 2 + pj) 



where P is the projector onto the 2x2 upper left corner of the N c x N c matrix of 
the instanton orientation. Just as in the case of the zero modes the near mass shell 
asymptotics of these propagators coincides with the perturbative propagators up to 
factors 1 / vTI: 
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Sj{x,z) 



1 -z 



Si(z,x) 



1 z 
/TL.2it 2 z 4 



(10) 



Thus amputation of the external legs in the near mass shell limit, denoted by a cross 
in the diagrams in Fig.l, leaves the spinors k and (p for the zero mode legs, see eq.(||), 
and a factor 1/vTT for the propagator legs from (|T0|). The leading term in the operator 
product expansion of the current product in the instanton background is then given by 



Jd^e^{j,(x)Uy))\^ T f = /dAe" A J % d(p) dx dU (An 2 p 3 )"^ 1 JJ [(x i¥ >)(/#i)] 



ijtu,d 



xl (An p ) tr {cr lx Si[x,y)a v Si{y,x)) + tr (a IM S I (x, y)a u S T (y, x)) (x«^)(«'0«)(Xd¥')(^d) 



2np 



3/2 



(Xu ( p){^Po(y)cruSi(y,x)a^ u )(xdf){^d) 



(11) 



+ (Xu^p)(Ktpu)(Xd^nSi(x,y)a v Ko(y))(Kip d ) 



+ (u^d,p<^v,x^y) 



+ 



{ (Xu<P>) (VO (x)o^u) (Xdf) (fO (y)<Tv1pd) + [XuGvKO (y) ) (^u)(Xd^^ (x) ) (Klp d ) } 



+^{{XuVvKo{y))(Vo(x)(T f ,'ll) u )(Xu<p)(K1pu) + (Xd^^O (#)) (<f0 (y)Vv4>d) (XdV) («^d) } 



Only the last term coming from the last diagram of Fig.l differs by its sign for the 
vector and axial vector currents. The integration over the instanton size is performed 
with the instanton density 



d(p) 



(N c -l)\{N c -2)\ 



-N c C2+n,fC3 



27T 

a(p) 



2N C 



-2n/a{p) 



(12) 



where c\ = 0.466 and the constants 02,03 take the values C2 = 1.54, C3 = 0.153 in the 
MS scheme. Now it is only a matter of patience to integrate over the instanton position 
and size and to take the Fourier transform. A typical p-integral has the structure 



dp 



,25+1 ln E 



p 2 A 2 



Ki(pq) 



(13) 



where K\ is a modified Bessel function and B, E are some numbers. The logarithmic 
size- dependence stems from three sources: from the two- loop running of a(p) in the 
exponent of the instanton density, from the preexponential factor a(p)~ 2Nc and from the 
anomalous dimension of the operators in the effective Lagrangian which are normalized 
at the scale fi = 1/ p.[] To perform this integral, we first calculate a similar integral 

t For consistency, we adjust the values of the one-loop and the two-loop A to reproduce the same 
value of a{m T ), which is the only parameter appearing in the final formulas. 
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with the logarithmic factor substituted by a power-like one (pA) 2e , and introduce the 
projection operator 

where the contour C wraps around the negative real axis. For example, for integer E, 
or general E and lnx > 0, one has V% {x e } = \n E x. For arbitrary real E this operation 
produces an asymptotic series 

Mffi'^ffis^^ (15> 

with Z(e) = z n£ n being the result of the p-integration. For natural E the series 
terminates and yields an exact answer for integrals of the type (|T3|). For general real 
E the asymptotic expansion in flI5| ) yields a good approximation provided ln(g 2 /A 2 ) is 
sufficiently large, and provided the integral is dominated by p < A -1 . 

Using the projection operator defined in (0), we can write down the instanton con- 
tribution to the two-point correlation function of (axial) vector currents in the back- 
ground of vacuum fermion fields in the compact form 



c e -2,/a( q) J dU Y[ (x^)(^)(-/3oa(g)) s 



s q 2 J \a(q) 



2 



+ (X«¥ 3 )(«[^(U - q^u)(Xd<p)(K[qa u - q u ]^d) 
+ {Xu[o v q - qu}^){Rip u ){Xd[cr^q ~ q^}f){^d) 
+ A (XuWuCL - qv]v)(Ki) u )(xuV)(K[<l d ti - qMu) 
+ A (xdWf.q - g A J^)(«VO(Xd | P)(«[?^ - qiM>d) 



(16) 



Here we defined 



B= b -±^, b = -4n(3 , E = 2N c - 2 -^ + 7 , (17) 

where (3 , (5\ are the first two coefficients of the /3-function, /3(a) = /i 2 <9/<9p 2 a = 
(3 n ot n+2 . The quantity 7 stands schematically for the anomalous dimensions of the 
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six-quark operators - that should be normalized at q - appearing in the square brackets 
in eq . (|16D . An overall factor of two, which accounts for an equal contribution in the 
background of an antiinstanton, is included in (p!6|). 

A comment is in order, concerning the integration over the instanton size. The 
integration over p contains a divergent term ~ {S^uq 2 — q^qv) V? 4 / dpp 2B ~ 5 . This term 
comes from large distances and has to be identified with an instanton contribution to 
the matrix element of the operator aG^G^. It will be neglected in the following. 
The remaining terms are convergent and represent the instanton contribution to the 
coefficient functions of the 2nj-quark operators shown in fllEp. The distinction between 
contributions to the matrix elements and to the coefficient functions is unambiguous. 
The relevant terms can easily be distinguished by their g 2 -dependence. 

The answer in (|T6|) has to be averaged over the physical vacuum. For numerical 
estimation, we assume factorization of the vacuum expectation values of the relevant 
six-quark operators, which allows to express them in terms of the quark condensate (qq). 
Upon factorization, the second and the third term in the square brackets in ( |16|) vanish. 
Since we are working in the chiral limit m u = rrid = m s = 0, we use the S , f/(3)-flavour 
symmetry (uu) = (dd) = (ss), and take the numerical value (qq) = — (240 MeV) 3 at 
the scale 1 GeV. Extracting the Lorentz structure, n^(g) = (q^q u — S^q 2 ) U ud (q 2 ), and 
taking the trivial integration over the instanton orientation, we arrive at 



2 \ "/ 



UJ W (N c -l)\(N c -2)\ 



-N c C2+rtfC 3 



% 1 r 



n 

i=l 



(liQi)M 



-a 



-2n/a(q) 



1 



a(q) 



27V c + 7 



x {-P a(q)) E V% \ e- 6 ^ 001 ^ F(e) 



1 

2B - 1 - 2e 



(19) 



For the anomalous dimensions of the six-quark operators factorization amounts to tak- 
ing 7 = Anf/b, where 4/b is the anomalous dimension of the quark condensate. However, 
since factorization is known not to be consistent with the renormalization group, we 
should rather allow 7 to vary in a certain range, ascribing the variation of U.f^ st to the 
uncertainty inherent to the factorization hypothesis.^] 

Corrections to eq.(^) come from (i) higher-dimensional operators, suppressed by 
powers of (pA) 2 and (ii) exchange of hard particles, suppressed by powers of a (/>*), 
where p* is the average size of instantons that contribute to the correlation functions. 
From now on we shall discuss the case rif = N c = 3,q = m T , which is the case of physical 
interest in r decays, and regard n^ si as a function of a(m T ). Then b = B = 9 and 

* Strictly speaking, in presence of instantons, the renormalization group equations for (2n/)-quark 
operators are more complicated and include mixing with the unity operator. This can be checked by a 
direct calculation of the (anti)instanton contribution to the vacuum expectation value ({qq) nf ), which 
proves to contain a logarithm of the renormalization scale. Taking into account this mixing corresponds 
to the calculation of correlation functions in the instanton-antiinstanton background in the spirit of 
fl4| , but does not seem appropriate for our present purposes since the size of the antiinstanton in this 
calculation turns out to be unacceptably large. 
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E = 34/9. To estimate the value of a(m T ), at which the above mentioned corrections 
become of the same order as the contribution kept in (|19D , we should first find p*. The 
distribution of instanton sizes relevant to the term proportional to A in ( |T9"D is given by 
the integral 

^ 2B+I ^ /*• (4=* (^) fell w 



P J p*A." J \y/UU \\fUUj UU \yJUU 



with u = 1 — u. The w-integration is sharply peaked at u = 1/2, which allows to put 
it = 1/2 for examination of the p-distribution. We plot the real part of this distribution 
in Fig. 2 as a function of a(m T ) and use a(m T ) = (—P \nml/A 2 )~ 1 to relate A and 
a(m T ). 

The distribution clearly shows two peaks. The first one is located at values p ~ 
(4 — 5)/m T , the precise coefficient depending weakly on E and a(m T ), and corresponds 
to the effective value of the instanton size p* ~ 400 MeV. The second one corresponds 
to the integration region p > 1/A, and its contribution should be small enough (say, 
by factor 5) in order that the calculation makes sense. Combined with the requirement 
that the value of the coupling at the effective scale is not too large, say a(p*) < 1, this 
restriction yields the critical value of a(m T ), above which the instanton contribution is 
ill-defined 

a(m T ) cr w 0.32 , (21) 

even if Hf^ st is still small. This value of the coupling lies well within the interval under 
discussion a{m T ) ~ 0.28 — 0.38 Note that since p*m T ~ 4 — 5 is a large number, 
the effect of including the logarithms lnl/(p 2 A 2 ) into the integration is important, and 
yields a suppression by roughly a factor 

(Mi/^) £ „ (*JpA) E „ 10 - (22) 

^In(m2/A 2 ); UW/ 

for typical values of E and a(m T ). Technically, this effect comes from summation of 
higher-order terms in the expansion in (|15D . 

3. We are now in the position to calculate the instanton contribution to the r decay 
widths. As usual, they are normalized as 

_ T(t~ -> v T + hadrons) 
T{t -> v T e v e ) 

Decomposing the correlation functions as 

^u,v/a(i) = (9m?" ~ V? 2 ) n (i),y/A(^ 2 ) + n (o),VM(? 2 ) ( 24 ) 
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and using analyticity in the cut q ,2 -plane, one can express the total hadronic width R T 



as 

2 



ds f s 

R T = 67TZ / ;r- 1 



m 2 V mi 



1 + 2 ^2 ) %)(«) + n (0)( s ) 



(25) 



where s = —q 2 and 

n (J) (s) = |K d | 2 (n$,vW + n (jW*)) + IK,| 2 (nfj )iV («) + n?j)») (26) 

with the CKM matrix elements V u d,V us . Experimentally the contribution from the 
strange current, R Tj s, can be separated according to the net strangeness of the final 
state. The nonstrange contribution can be further resolved into vector and axial vector 
pieces, R T y an d R t .a, according to whether the final state contains an even or odd 
number of pions. Thus, one defines 

Rr = Rr,v + Rt,a + R T ,s = 3 (\V ud \ 2 + \V US \ 2 ) {1+5), (27) 

R T ,v/A = l\V ud \ 2 (l + 5 v/A ), (28) 

separating the "Born" term from the QCD corrections 5 (a small multiplicative elec- 
troweak correction is understood), that are of the order of 20%. We will now determine 
the instanton contribution S inst ^ Syf^ to the QCD corrections. In the massless limit the 
correlation functions are transverse, see eq. (|T9|) , i.e. il^ s *(g 2 ) = 0. The contour integral 
in Q25D is then taken including logarithms of s/A 2 exactly, using the same trick as above 
for the integration over the instanton size. This produces immediately 



xinst 
°V/A 



-7T • 4 



B 



2n 2 



(2vr) 6 ci e - 3c2+3c3 J] 



i=l 



-a 



m-i 



(/^) 7/3 



-2ir I a(m T ) 



(29) 



(— ) (-(3 a(m T )) E V% \ e-^d)) F{e) H{e) 



2B - 1 - 2e 



± 1 



which is our final result. The upper sign holds for vector currents, the lower sign for 
axial vector currents. The effect of the contour integration - up to constant factors - 
is contained in the function H(e), 



H(e) = sin tt(B - e) 



1 



1-B+e 3-B 



+ 



A-B + e 



(30) 



Note that the expansion of H(e) starts with order e, since B is integer. Therefore the 
transition from the correlation functions to the decay widths suppresses the instanton 



8 



contribution by one power of a(m T ). With (ss) = (uu) the instanton contribution to 
the strange decays equals that to the nonstrange decays. Then 

§ inst = 1 (gvnst + ginst^ _ i_ gmst _ ^ 

The term differing in sign for the vector and axial vector widths in (^) cancels in the 
total width, leaving only the first term with the small coefficient ~ 1/(2-8 — 1). Thus 
there is a strong cancellation in the sum of vector and axial vector contributions, which 
is similar in effect to the cancellation in the contribution from four-quark operators 
in the standard SVZ expansion [|I|]. However, this cancellation relies strongly on the 
factorization approximation, which eliminates the second and third term in (|T6|) having 
much larger coefficients. For this reason, the instanton contribution to the total width 
is much more sensitive to deviations from factorization and therefore less reliable nu- 
merically. Since 5y st = —8™ st + 2S mst , the axial vector contribution is essentially of the 
same magnitude as for the vector current but with an opposite sign. In Fig. 3 we plot 
5y st and 6 mst as a function of a(m T ) for three choices of the anomalous dimension of 
the (qq) 3 operator 7 = 2/3, 4/3, and 6/3, corresponding to E =28/9, 34/9, and 40/9, 
respectively. (The middle value corresponds to factorization.) At a(m T ) = 0.32 we get 

§ mst _ _ 6 inst = Qm _ Qm ^ S inst _ Q.002 _ 0.003 . (32) 

At larger values of the coupling a(m T ) ~ 0.34 — 0.36 the instanton contribution blows 
up, but our calculation is no longer under control. Thus the instanton contribution is 
essentially of the same size as the contribution from dimension-6 operators ~ (qq) 2 /q 6 

[0, 

<S£= 6 = 0.024 ± 0.013 , <5^ =6 = - (0.038 ± 0.020) , 5 D=6 = - (0.007 ± 0.004) , (33) 



which are the largest power corrections in the standard approach. 

4. We performed an explicit calculation of the instanton contribution to the r-decay 
widths, which corresponds to the exponential correction to the coefficient function in 
front of the six-quark operators in the operator product expansion of the relevant cor- 
relation functions, and which is distinguished by the chiral properties of the instanton- 
induced effective vertex. Contrary to previous claims, we find a sizable instanton con- 
tribution, which for a(m T ) w 0.32 is of the same order of magnitude as standard 
nonperturbative corrections induced by nonzero vacuum expectation values of the op- 
erators of lowest dimension. An inherent uncertainty of our calculation is about a factor 
of two for vector and axial vector channels, and maybe larger for the total width owing 
to strong cancellations. For still higher values of a(m T ) we cannot make quantitative 
statements, but there is no reason to expect that the instanton contribution becomes 
small. 



9 



Unfortunately, as far as numerical values are concerned, our result is of limited 
practical importance, since the effective instanton size has proved to be large, ~ 
400 MeV. At such low scales, the instanton density is modified strongly by external 
gluon fields. The effect is conceptually quite similar to the effect of the quark condensate 
in ([[]), and the lowest-order correction reads JTTJ] 



d eff(p) = d(p) 



1 + sJ{p) (tt G2 



(34) 



where (a/nG 2 ) ~ 0.012 GeV 4 is the gluon condensate. With the effective instanton 
size 400 MeV, the correction term in brackets in (^) is several times larger than unity, 
indicating that the effect is important. For the case of r-decay, this means that instanton 
contributions to the coefficient functions of operators of higher dimension, including 
gluon fields in addition to six quark fields, are likely to be larger than the one considered 
in this paper, so that the expansion in powers of l/m T fails. One may speculate that 
this expansion is at best an asymptotical one, and the contributions ~ 1/m* 8 under 
discussion are already in the region where the series starts to diverge. Anyhow, in 
agreement with the old wisdom [jllj, we find that instanton calculations are hardly 



useful for quantitative estimates of nonperturbative effects, but rather can be used to 
indicate a scale at which the power expansion breaks down. 

Thus, our conclusion is partly pessimistic. On the one hand, we find it difficult to 
justify the use of the operator product expansion in finite energy sum rules at the scale 
of the r mass, since contributions of higher orders are of the same order as leading power 
corrections. Our calculation essentially supports the old philosophy of SVZ ||, who 
introduced Borel sum rules in which higher-order power corrections are suppressed by 
factorials. In our case Borel transformation would introduce a factor 1/8! ~ 2.5 • 10~ 5 
and render the instanton contribution completely negligible. In the practice of the 
numerous QCD sum rule calculations this old argumentation has partly been forgotten. 
The reason is that higher-order corrections have never been calculated explicitly, and 
as far as leading-order power corrections are concerned, Borel and finite energy sum 
rules give very similar results, within the typical 10% accuracy. It is the requirement of 
a very high precision, which makes mandatory to reconsider the theoretical accuracy of 
the sum rule program for r-decays. In the light of our result a simultaneous extraction 
of a(m T ) and nonperturbative parameters from r decays J| is hard to justify. 

On the other hand, all nonperturbative contributions to the total r-decay width 
prove to be small, less or of order 1%. Although the leading power corrections cannot 
be guaranteed to improve the accuracy, they can well indicate an intrinsic uncertainty 
of the theoretical prediction. Thus, it is plausible to expect an accuracy of the per- 
turbative prediction for the r hadronic width of order 1%, but the situation could be 
substantially worse for the exclusive vector and axial vector channels. This allows for a 
determination of the QCD coupling within 10% accuracy, which is competitive to the 
current accuracy of the determination of a at the scale of the Z-boson mass, but, in 
difference to the latter, cannot be substantially improved. 
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Figure Captions 



Fig.l The four types of diagrams that contribute to the Wilson coefficient of six-quark 
operators. The instanton is depicted as a 2n/-quark vertex with only the flavours 
u and d shown explicitly. A circled cross denotes the insertion of the current. 

Fig. 2 Distribution of instanton sizes as a function of a(m T ). The vertical scale is 
arbitrary. 

Fig. 3 Instanton contribution to the vector decay channel (a) and the total width (b) 
for three values of "anomalous dimensions", E = 28/9,34/9,40/9, as a function 
of a(m T ). 
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